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Abstract 

It is shown that in the Love-KirchhofF plate theory, an edge wave 
can travel in a circular thin disk made of an isotropic elastic ma- 
terial. This disk edge wave turns out to be faster than the classic 
flexural acoustic wave in a straight-edged, semi-infinite, thin plate, a 
wave which it mimics when the curvature radius becomes very large 
compared to the wavelength. 

1 INTRODUCTION 

When an elastic half-space is subjected to a point excitation, it responds 
by generating waves which become surface (Rayleigh) waves in the far-field. 
When a cymbal is hit on its edge, it produces a long-lasting sound: why not 
conjecture that this phenomenon is due to the propagation of an acoustic 
wave localized near the edge of the cymbal? 

The study of vibrations of thin circular plates with free edges dates back 
to at least 1910 with a paper by Airey [T], and is still a very active area of 
research (see for instance Touze, Thomas, and Chaigne [2]). In contrast there 
seems to be almost no contribution on the subject of waves in these plates, 
whereas there exist a wealth of papers on waves over the surface of cylinders 
with infinite extent. That scarcity might be explained as follows. In the 
latter case (cylinder), the asymptotic behavior of a wave when the radius of 
the cylinder becomes very large (in comparison to the wavelength) is well- 
known: it is that of the celebrated Rayleigh wave in a half-space [3J. In the 
former case (disk), the asymptotic behavior is not as well-known: it is that of 
a wave propagating along the straight, free edge of a thin Kirchhoff plate [1]. 
Hence Cerv ^ studied waves traveling on the surface of a thin disk, but used 
plane-stress equations, where the thickness of the plate does not appear: as 
a result the asymptotic behavior of his waves is that of Rayleigh waves. 
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This paper shows that the appropriate resolution of the classical small- 
deflection equations of motion in a circular disk causes no particular problem. 
It is conducted with a view to interpret some patterns of motion observed 
near the edge of a stricken cymbal [6] (other experimental results are found 
in [7]). 

Other related problems are left open but should pose no particular diffi- 
culty to solve; they include the following, (i) The complimentary problem of 
acoustic waves propagating along the edge of a circular hole in a thin plate, 
with obvious possible applications in non-destructive evaluation, see Figure[Il 

(ii) A circular disk made of one material in rigid or slippery contact with a 
holed plate made of another material; there the asymptotic behavior is that 
of contact fiexural waves in thin plates [S] and the applications are also in the 
realm of non-destructive evaluation (of the strength of rivets for example). 

(iii) Waves traveling in a thin fiat annulus; there the asymptotic behavior is 
that of fiexural waves on a narrow plate p] and an exotic application would 
be the transportation of planar heavy objects by acoustic levitation [lOj on 
circular transporting systems, (iv) Another possibility resides in linking the 
present results to the experimental work of Shaw [jjj on Barium Titanate 
disks. 

2 BASIC EQUATIONS 

The equations of motion for a Love-KirchhofF plate are [T2] 

{VshyV^w + Q{l-u)w = 0, (1) 

where w is the deflection, h is the thickness of the plate, u its Poisson ratio, 
and Vs its shear wave speed. Now consider an acoustic fiexural wave traveling 
along the circumference of a circle (axis normal to the plate, radius R) in the 
form 

w = W{r)cos{kRe -ut), (2) 

where r, 6 are the in-plane cylindrical coordinates, k is the wavenumber, and 
uj the frequency. Introducing the dimensionless "speed" a, defined by [13] 

the equations of motion can be put in the form 

B+D^W = D_D+W = 0, (4) 
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Fig. 1: Wave localized near the edge of a disk or a hole. 



where D± are the following differential operators: 

D± := dVdr^ + (1 /r)d/dr - k'^R^/r'^ ± afcl (5) 

Clearly, the general solution to this differential equation is a linear com- 
bination of Bessel functions; specifically, 

W{r) = 7iJp(A/aA;r) + 'y2lp{y/akr) + 'ysYp{y/akr) + ^4Kp{y/akr), p := kR. 

(6) 

Here, the ji are constants and the labeling of Bessel functions and modified 
Bessel functions follows the conventions of Abramowitz and Stegun [Tl]. In 
particular, 

D+Jp(v^A;r) = D+Yp(v^A;r) = DAp{^/akr) = D_Kp{^/^kr) = 0. (7) 

Then for the edge of the circle to be free, two boundary conditions must 
be satisfied [12]. First, the edge shear must vanish, so that 

d'^w ^ /Idw ^ 1 d'^w\ t i? (8) 

Second, the rate of change of the twisting moment must also vanish, so that 

9 {1 — ly) d d^w 1 dw\ O t R (9) 

dr r 06 \r drdO 89 J 

Finally, the displacement at 9 is the displacement at ^ + 27r, 9 + 47r, etc., 
and so by (I2l), p = kR must be an integer. In other words, the wave can 
only travel under the condition that its wavelength is an integer fraction of 
the circle perimeter 2tiR. 

3 CIRCULAR DISK 

For a wave on the edge of a circular disk, take 73 = 74 = in (I6l) to avoid 
singularities at the center of the disk. Then the boundary conditions (IHl) and 
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© reduce, after much manipulation and use of ([7]), to 



^piVap) [(1 - u)jp{a)/p - (1 - - a)] 71 

+ IpiV^p) [(1 - 'y)tpia)/p + (1 - z/ + a)] 72 = 0, (10) 



and 

^piVap) [{l-iy + a)jp{a) - (1 - iy)/p] 71 

+ IpiV^p) [{l-i^- a)ip{a) - (1 - u)/p] 72 = 0, (11) 

respectively. Here the functions jp{a) and ip{a) are defined by 



jp(a) := -^/a /. _ , ~ y/l — a as p — * 00, 
Ip(V"p) 



ip(a) := ^/a J' ~ vT+a as p ^ 00, (12) 

IpWap) 

and their asymptotic behavior as p becomes large is found from Debye series 
expansions [H]. For non-trivial solutions to exist, the following dispersion 
equation must be satisfied: 

{1 — u — a)'^ip{a) — (1 — z/ + a^jp^a) 

+ 2a[l + tpia)jp{a)]il - u)/p - [tpia) - jp(a)](l - uf/p^ = 0. (13) 

Not surprisingly, the dispersion equation for disk edge waves is the same as 
the frequency equation [H [16] derived for the vibrations of a circular disk 
with free edges, simply because the governing equations and the boundary 
conditions happen to coincide in the present context for waves such as ([2]) 
and for vibrations in the form w = W{r) cos{kR6) cos{ujt), say. The inter- 
pretation is however different. In particular, when the radius of the disk 
becomes large with respect to the wavelength, p 00 and this equation be- 
comes that of the edge wave on a semi-infinite plate discovered by Konenkov 
|1] and hereafter called the "straight-edge wave" , 



u — a 



Vl + a = (l-z/ + a)Vl-a. (14) 



Note that in general, the dispersion equation (fT3l) gives several roots, 
numbered by an integer parameter s, say. The smallest root (s = 0) corre- 
sponds to a solution which is never zero over the surface of the disk; the next 
root (s = 1) corresponds to a solution which is zero once over the surface of 
the disk, and so on. Hence s denotes the number of nodal circles [16\ . 
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Table 1 displays the values of ot for various Poisson ratios and integer 
values of p at s = 0; clearly, a tends to the straight-edge wave speed as 
p oo (last column, see also Thurston and McKenna [IH]); also, the disk 
edge wave is always faster than the straight-edge wave (note that similarly, 
the wave propagating on the free surface of a cylinder is faster than the 
Rayleigh wave propagating over the free surface of an elastic half-space |15].) 
For the one- mode {p = 1) wave, a is nearly twice its value for the straight- 
edge wave {p = oo) and then it decreases rapidly as p grows. Multiplying the 
values displayed by p"^ gives the quantity computed and collected by Leissa 
|16j : the reader is also referred to that book for values when s > 1. 



Table 1. Values of a for various v and p (s = 0) 



\ p 


5 


10 


20 


50 


100 


oo 


0.1 


1.43125 


1.30545 


1.20294 


1.11354 


1.07218 


0.999985 


0.17 


1.40355 


1.28864 


1.19284 


1.10814 


1.06862 


0.999852 


0.2 


1.39035 


1.28028 


1.18757 


1.10510 


1.06656 


0.999697 


0.25 


1.36639 


1.26460 


1.17736 


1.09907 


1.06227 


0.999177 


0.3 


1.33980 


1.24654 


1.16513 


1.09144 


1.05668 


0.998102 


1/3 


1.32048 


1.23302 


1.15570 


1.08533 


1.05209 


0.996899 


0.4 


1.27763 


1.20198 


1.13326 


1.07013 


1.04031 


0.992639 



For the purpose of comparison, recall that Viktorov [IS] studied waves 
over the convex surface of a full cylinder, plotted their variation with depth 
for a cylinder made of brass, at p = 5,41,cxd, and found that they 'decay 
somewhat more rapidly with distance from the surface than for a Rayleigh 
wave, the rate of decay being faster for smaller p = kR\ Figure 2 shows 
that the depth profiles (normalized with respect to the wavelength) at p = 
5, 41, 200, oo and s = for a brass disk follow a similar trend. Note however 
that the edge waves on a disk decay much more rapidly than the straight-edge 
wave, which is known to be weakly inhomogeneous [17] . 

Schedin, Gren, and Rossing [6] used double-pulsed television holography 
to record the response of a cymbal made of brass (z/ = 0.34) when it is struck 
near its edge. They noted that 'high-frequency waves are apparent along the 
edge' and these can be clearly seen on their pictures. A rough estimate gives 
a number of about 12 wavelengths in a ttR/S sector, indicating that p is about 
200. Figure 3 displays the dispersion curves for brass for s = 0, 1, 2, 3, up to 
p = 200. Here, the discrete nature of the plots has been over-ridden for ease 
of reading. Finally, the profiles of the disk are shown on Figure 4 at p = 200 
for s = 0,1,2; there the distance was normalized with respect to the radius 
of the disk and the elevation at the edge of the disk was taken as half the 
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Fig. 2: Depth profiles for a tliin disk made of brass (z/ = 0.34). 

radius (exaggerated scale). For s = (no nodal circle), the amphtude decays 
rapidly away from the edge, and is almost zero at 0.9-R. For s = 1,2, the 
amplitude of the wave reaches a maximum which is greater than its value at 
the edge (amplification); it is also less localized than for s = 0, being almost 
zero at about 0.8R. Seen from above, the case s = 2 (two nodal circles) will 
display one shaded region and one bright region close to the edge, a pattern 
which is somewhat similar to that recorded by Schedin, Gren, and Rossing 
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Fig. 3: Dispersion curves for a disk made of brass (z/ = 0.34). The dashed 
hne represents the speed of the straight-edge wave. 



0.2 n 
0.1 



s = 



s = 2 



r/R 



0.2 



0.4 



0.6 






-0.1 
-0.2 

Fig. 4: Profile of a disk made of brass [v = 0.34) when p = 200, s = 0, 1, 2. 
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